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We present a search for B decays to charmless final states involving charged or neutral ao mesons. 
The data sample corresponds to 89 million BB pairs collected with the BABAR detector operating 
at the PEP-II asymmetric-energy B Factory at SLAC. We find no significant signals and determine 
the following_90% C.L. upper limits: B(B° -> a ( 77r+) < 5.1 x 1(T 6 , B(B° -» a^ K+) < 2.1 x 1CT 6 , 
B(B~ -> ao K°) < 3.9 x 1(T 6 , B{B+ -» a^ir+) < 5.8 x 1CT 6 , B(B+ -» a° K+) < 2.5 x 1(T 6 , and 
B{B° -» ag-K" ) < 7.8 x 1(T 6 , where in all cases indicates the product of branching fractions for 
B — > aoX and a rjn, where X indicates 7f or 7T. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



We report results on measurements of B-meson decays 
to charmless final states with ao(980) mesons 0. Both 
experimentally and theoretically, most work in charm- 
less two-body B decays has involved states with only 
pseudoscalar and vector mesons. The only charmless B 
decay involving scalar mesons that has been observed 
is B — > f (980)K [2(. There have been no previously 
published searches for B decays to final states with ao 
mesons. In this paper we search for the decays B — > ao7r, 
B — > a$K , and B — > clqK for both charged and neutral 
ao mesons. These measurements should provide informa- 
tion both for B decays to scalar mesons and the nature 
of those mesons. 

Some specific predictions can be made for the decays 
B — > fto^ if factorization is assumed and if the decay 
is a tree or penguin (loop) process. The dominant such 
process is shown in Fig. ^a). The companion tree pro- 
cess, shown in Fig. QJb), is expected to be greatly sup- 
pressed, since the virtual W cannot produce an ao meson 
[4J. This is a firm prediction of the Standard Model be- 
cause the weak current has a G-parity even vector part 
and a G-parity odd axial- vector part. The latter can pro- 
duce an axial-vector or pseudoscalar particle while the 
former produces a vector particle, but neither can pro- 
duce a G-parity odd scalar meson. Penguin processes 
such as shown in Fig. ^c) are allowed, but are sup- 
pressed relative to the tree processes. Thus the decay 
B — ► clq7t^~ is expected to be "self-tagging" (the charge of 
the pion identifies the B flavor). The decays with a kaon 
in the final state should be dominated by penguin pro- 
cesses (Fig. |3i); however, there is a cancellation between 
two terms in the penguin amplitudes for these decays || , 
which leads to a prediction that the branching fraction 
should be rather small. The diagrams for neutral B de- 
cays involving a§ mesons are similar to those shown in 

pig. ni 

The nature of the ao is still not well understood. It 
is thought to be a qq state with a possible admixture 
of a KK bound-state component due to the proximity 
to the KK threshold The ao mass is known to be 

about 985 MeV and the dominant decay mode is ao — > t/tt 
0, which is the mode used in the present analysis. A 
recent analysis || that uses this ijir decay channel finds 
a Brcit-Wigner width of (71 ± 7) MeV, with no better 
fit obtained when the more correct Flatte shape || is 
used. Also since the branching fraction for ao — * rjir is 



B" 




d d d d 



(c) (d) 

FIG. 1: Feynman diagrams for decays involving charged 
ao mesons: (a) dominant and (b) G-parity-suppressed tree 
diagrams for B° —* a^^, (c) penguin diagram for the 
same decay mode, and (d) penguin diagram for the decay 
B°^aWK + . 



not well known, we report the product branching fraction 
B(B — > a A") x B(a Q — > rjir), where X indicates K or ir. 

The results presented here are based on data collected 
with the BABAR detector ^ at tne PEP-II asymmetric 
e + e~ collider located at the Stanford Linear Accelerator 
Center. An integrated luminosity of 81.9 fb _1 , corre- 
sponding to 88.9 ± 1.0 million BB pairs, was recorded 
at the T(4S) resonance (center-of-mass energy y/s — 
10.58 GeV). 

The track parameters of charged particles are mea- 
sured by a combination of a silicon vertex tracker, with 
five layers of double-sided silicon sensors, and a 40-layer 
central drift chamber, both operating in the 1.5-T mag- 
netic field of a superconducting solenoid. We identify 
photons and electrons using a CsI(Tl) electromagnetic 
calorimeter (EMC). Further charged particle identifica- 
tion (PID) is provided by the average energy loss (dE / dx) 
in the tracking devices and by an internally-reflecting, 
ring-imaging Cherenkov detector (DIRC) covering the 
central region. 

We select ao candidates from the decay channel ao — ► 
i]tt with the decays r/ — > 77 (?7 77 ) and 77 — ► 7r + 7r~7r° 
(77371- )■ We apply the following requirements on the in- 
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variant masses (in MeV) relevant here: 500 < m 77 < 585 
for ?y 77 , 535 < < 560 for 773^ , 120 < m 77 < 150 

for 7T°, and 775 < to,™ < 1175 for a — > t]tt. These 
requirements are typically quite loose compared with 
typical resolutions in order to achieve high efficiency 
and retain sufficient sidebands to characterize the back- 
ground for subsequent fitting. We reconstruct K® can- 
didates through the Kg — > 7r + 7r~ decay; to obtain a 
low-background, well-understood K® sample, we require 
488 < m^Tr < 508 MeV, the three-dimensional flight dis- 
tance from the event primary vertex to be greater than 
2 mm, and the angle between flight and momentum vec- 
tors, in the plane perpendicular to the beam direction, 
to be less than 40 mrad. 

We make several PID requirements to ensure the iden- 
tity of the pions and kaons. Secondary tracks in rj^ can- 
didates must have measured DIRC, dE/dx, and EMC 
outputs consistent with pions. For the decays B — ► aoh + 
0, where h + indicates a charged pion or kaon, the par- 
ticle h + must have an associated DIRC signal with a 
Cherenkov angle within 3.5 standard deviations of the 
expected value for either a ^ or K hypothesis (we de- 
scribe below the separation between the two hypotheses). 

A _B-meson candidate is characterized kinemat- 
ically by the energy-substituted mass ?ties 
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Po • Pb) 2 /Eq — p 2 B ]i and energy difference 
AE = E* B — \\fs, where (Eb,Pb) and (i?o,Po) are 
the four vectors of the B-candidate and the initial 
electron-positron system, respectively. The asterisk 
denotes the T(4S) frame, and s is the square of the 
invariant mass of the electron-positron system. The 
AE (toes) resolution is about 40 MeV (3.0 MeV). We 
require \AE\ < 0.2 GeV and 5.2 < m ES < 5.29 GeV. 

Backgrounds arise primarily from random combina- 
tions in continuum e + e~ — > qq (q = u, d, s, c) events. We 
reduce these by using the angle 9t between the thrust 
axis of the B candidate in the T(45 l ) frame and that 
of the rest of the charged tracks and neutral clusters in 
the event. The distribution of | cos#t| is sharply peaked 
near 1.0 for combinations drawn from jet-like qq pairs, 
and nearly uniform for i?-meson decays. We require 
|cos#t| < 0.9 for the aoK® decay modes. Based on a 
Monte Carlo study in which the relative branching frac- 
tion uncertainty is minimized, we tighten this require- 
ment for the higher-background a^h channels: 0.8 for 
o-o (V3Tr)h + , 0.7 for Oq (r; 77 )/i + and aQ(r]s n )h + , and 0.6 
for Oq(?7 77 )/i + . We also use, in the fit described below, a 
Fisher discriminant T that combines the angles with re- 
spect to the beam axis of the B momentum and B thrust 
axis (in the T(4S) frame), and moments describing the 
energy flow about the B thrust axis (ll| |. 

For the 77 — > 77 modes we use additional event- 
selection criteria to further reduce backgrounds from 
charmless B decay modes such as B — > K* 7 and B — > 
rjK* . We require |cos#j ec | < 0.86, where 6J CC is the rj 
decay angle, the angle of the photons in the 77 rest frame 



with respect to the boost direction from the B to that 
frame. We also require cos6^° c < 0.8, where #^° c is the 
ao decay angle, defined similarly to 0^ ec , with sign such 
that high-momentum rj mesons populate the region near 
+1. These additional requirements reduce the BB back- 
ground by a factor of 2-4, depending on the decay mode. 
From Monte Carlo (MC) simulation [12J we estimate that 
the residual charmless BB background is less than one 
event for all decays except Oq (r] 71 )K° (the notation in- 
dicates the decay mode of the r\ used in reconstructing 
the ao) and etQ(?7 77 )ft, + , where we include in the fit a BB 
component, that we find to be less than 0.5% of the total 
sample in both cases. 

We obtain yields and branching fractions from ex- 
tended unbinncd maximum-likelihood fits, with input ob- 
servables AE, toeSj J 7 , m^Tti and for charged modes the 
PID variables S n and Sk] the last quantities are the 
number of standard deviations between the measured 
Cherenkov angle and the expectation for pions and kaons. 

For each event i, hypothesis j (signal, continuum back- 
ground, BB background), and, for the a$h + decays, fla- 
vor fc, we define the probability density function (PDF) 



V) k = VM^V^AEil, Sl\)V 3 (T l )V^). (1) 

The term in brackets for S pertains to the a$h + modes. 
The absence of correlations among observables (except 
between AE and S, which both depend on the momen- 
tum of the particle h + ) in the background Pj fe! is con- 
firmed in the (background-dominated) data samples en- 
tering the fit. For the signal component, we correct for 
effects due to the neglect of small correlations (more de- 
tails are provided in the systematics discussion below). 
The likelihood function is 



N 



£ = exp(-^F ife )n 



E 

3,k 



Y 3k V] k 



(2) 



where Yjk is the yield of events of hypothesis j and flavor 
k that we find by maximizing C, and N is the number of 
events in the sample. 

We determine the PDF parameters from simulation 
for the signal and BB background components, and ini- 
tial values of the continuum background parameters from 
(toes, AE) sideband data. We parameterize each of the 
functions P sig (m E s), V sig (AE k ), V 3 (T), andi^(5 fe ) 
with either a Gaussian function, the sum of two Gaus- 
sian functions or an asymmetric Gaussian function, as 
required to describe the distribution. The component of 
Vj{m ri T T ) which represents real mesons in the combi- 
natorial background is described with the same Brcit- 
Wigner parameters as arc used for signal. Slowly vary- 
ing distributions (clq candidate mass and AE for com- 
binatoric background) are represented by second order 
Chebyshev polynomials. The qq combinatoric back- 
ground in Toes is described by the function f(x) = 
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xy/ 1 — x 2 exp [— £(1 — a; 2 )] , with x = 2toes/-\/s and free 
parameter £; for Si? background, we add a Gaussian 
function to the quantity f(x). Large control samples of 
B — ► Z?7r decays of topology similar to the signal are 
used to verify the simulated resolutions in AE and toes- 
Where the control data samples reveal differences from 
MC, we shift or scale the resolution used in the likelihood 
fits. Examples of many of these PDF shapes from a very 
similar analysis are shown in Ref. [TT| . Additionally, 
the Breit-Wigner signal parameters for the do mass and 
width are determined from an inclusive dataset that is 
much larger than the sample used for this analysis. The 
widths are consistent with expectations from the natural- 
width values of Ref. || . 

In Table [I] we show for each decay mode the measured 
product branching fraction, together with the quantities 
entering into its determination. In order to account for 
the uncertainties in the background PDF descriptions, 
we include as free parameters in the fit, in addition to 
the signal and background yields, the principle param- 
eters describing the background PDFs: slopes for the 
polynomial shape for the AE and ao mass distributions, 
the parameter £ used in the Toes description, and three 
parameters describing the asymmetric Gaussian function 
for T . For calculation of branching fractions, we assume 
that the decay rates of the T(AS) to B + B~ and B°B° are 
equal We combine branching fraction results from 
the two 7] decay channels by adding the values of — 2 ln£, 
adjusted for a small fit bias (see below) and taking proper 
account of the correlated and uncorrelated systematic er- 
rors. 

In order to check the suitability of the PDFs for de- 
scribing the data, we show in Fig. the distribution 
of the likelihood ratio C(S)/[C(S) + C{B)\ for the full 
a o( 7 l-ry)h + sample, where C(S) and C(B) are the signal 
and background likelihood, respectively. Signal would 
appear near one in this plot but very little is seen because 
of the small signal yield. There is also good agreement 
for similar plots for the other samples. In order to show 
distributions of the main fit observables toes and AE, we 
require that this likelihood ratio be greater than a value 
that would optimize the branching fraction uncertainty, 
typically 0.9 for most samples. In Figs. 01 and 01 we show 
projections onto toes and AE of subsamples enriched 
with this requirement on the likelihood ratio (computed 
ignoring the PDF associated with the variable plotted). 

The statistical error on the signal yield is taken as the 
change in the central value when the quantity — 21n£ 
increases by one unit from its minimum value. The sig- 
nificance is taken as the square root of the difference be- 
tween the value of —2 In £ (with additive systematic un- 
certainties included) for zero signal and the value at the 
minimum, with other parameters free in both cases. The 
90% confidence level (C.L.) upper limit is taken to be 
the branching fraction below which lies 90% of the total 
of the likelihood integral (with systematic uncertainties 




L(S)/[L(S)+L(B)] 

FIG. 2: The likelihood ratio C(S)/[£(S) + C(B)] for 
a o{v~ii)h + ■ The points represent the on-resonance data, the 
solid histograms are from MC generated from background 
(dark shaded) and background plus signal (light shaded) 
PDFs. 




§.2 5.22 5.24 5.26 5.28 -$.2 -0.1 0.1 0.2 
m ES (GeV) AE(GeV) 



FIG. 3: Projections of the B-candidate m,ES and AE for (a, 
b) a^"/i + , and (c, d) agh + . Points with errors represent data, 
solid curves the full fit functions, dashed curves the back- 
ground functions (the peaking BB background component 
is negligible), and the dotted curve shows the kaon portion 
of the signal. These plots are made with a minimum require- 
ment on the likelihood and thus do not show all events in the 
data samples. 



included) in the positive branching fraction region. 

Most of the yield uncertainties arising from lack of 
knowledge of the PDFs have been included in the sta- 
tistical error since most background parameters are free 
in the fit. Varying the signal PDF parameters within 
their estimated uncertainties, we determine the uncer- 
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TABLE I: Signal yield, detection efficiency e, daughter branching fraction product dl^Oi significance (including additive 
systematic uncertainties, taken to be zero if corrected yield is negative), measured product branching fraction (see text), and 
the 90% C.L. upper limit on this branching fraction. 



Mode 


Yield 


e(%) Il^(%) Signif. 


S(10~ 6 ) 


UL(10" 6 ) 


a (r? 77 )7r 


18 -10 


18.8 


39.4 


1.3 


_ q-|-1.7 i r\ r\ 

z.6_ 1 5 ± u.y 




a ( r ?37r) 7r + 


io_ 8 


15.5 


22.6 


1.6 


o q+2.9 _i_ -1 n 




r.~ -r- + 








2.0 


1 3 ± U.7 


< 5.1 


a (?7 77 JK 


0+6 
*-4 


17.9 


39.4 


0.1 


U.U_q g ± U.o 




15+8 
1J -6 


14.9 


22.6 


1.1 


q I +2.5 ! I Q 












0.4 


0.4_ 8 ± 0.2 


< ^.1 


a (?? 77 )A 


- 12 -6 


21.4 


13.5 


0.0 


o 't+2.9 _i_ a n 

-3.7I 2 3 ±0.9 






n+ 7 

U -5 


15.8 


7.9 


0.5 


O 7+6.1 _i_ -i n 

z. i _ 4 4 ± i.y 




a K 








U.o 


— l.o_ 1 g ± U.o 


< 3.9 


ag(?7 77 )7r + 


17 +ii 

—9 


12.8 


39.4 


1.4 


± 1.2 




a o( 7 ?37r)7T + 




9.5 


22.6 


0.3 


i-2±|;l ± 1.7 




+ 

aoTT 








1.4 


2.&tx° 7 ± 1.0 


< 5.8 


a° ( VjJ )K+ 




12.4 


39.4 


0.3 


0.3± : i ± 0.4 








9.1 


22.6 


0.5 


1.9±|;| ± 2.5 




a° K+ 








0.4 


0.4+o 7 ± 0.3 


< 2.5 






15.0 


13.3 


0.5 


1.412.4 ± L2 








9.7 


7.8 


1.2 


6.6±£| ± 2.8 




a° K° 








1.0 


2.8+2 4 ± 11 


< 7.8 




m ES (GeV) AE(GeV) 

FIG. 4: Projections of the B-candidate rriES and AE for 
(a, b) Kg, and (c, d) ao-ffg. Points with errors represent 
data, solid curves the full fit functions, and dashed curves the 
background functions. These plots are made with a minimum 
requirement on the likelihood and thus do not show all events 
in the data samples. 

tainties in the signal PDFs to be 1-5 events, depending 
on the final state. The contribution to this uncertainty 
from the parameterization of the ao signal shape is small. 
We verify that the value of the likelihood of each fit is 
consistent with the expectation found from an ensemble 
of simulated experiments. 

Uncertainties in our knowledge of the efficiency, found 
from auxiliary studies, include 0.8%-iVt, 2.5%-iV 7 , and 
4% for a decay, where N t and N 7 are the number 



of signal tracks and photons, respectively. Our estimate 
of the number of produced BB events is uncertain by 
1.1%. The neglect of correlations among observables in 
the fit can cause a systematic bias; the correction for this 
bias (between —3 and +3 events) and assignment of the 
resulting systematic uncertainty (0.5-2 events) is deter- 
mined from simulated samples with varying background 
populations. Published data 6] provide the uncertainties 
in the _B-daughter product branching fractions (1-2%). 
Selection efficiency uncertainties are 0.5-3.5% for cos^t 
and 0.5% for P1D (for the a^h + modes). 

In conclusion, we do not find significant signals for 
these B-meson decays to states with ao mesons. The 
measured branching fractions and 90% C.L. upper limits 
are given in Table [I] Assuming r/ir to be the dominant 
ao decay mode, we rule out the predictions for the decay 
B~ — > %K° derived in Ref. Q. 
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